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Abstract. The quark-gluon plasma behaves as a relativistic viscous fluid in high-energy heavy 
ion collisions. I develop a causal dissipative hydrodynamic model at finite baryon density for 
RHIC and LHC to estimate the net baryon rapidity distribution. The net baryon number is 
found to be carried to forward rapidity by the flow, effectively enhancing the transparency of the 
collisions. This suggests that the energy available for the production of a hot medium could be 
larger than that naively implied from experimental data. Also the distribution can be sensitive 
to baryon dissipation as much as to shear and bulk viscosity. 



1. Introduction 

The quark-gluon plasma (QGP) [1] is a many-body system of quarks and gluons deconfined 
from hadrons at high energies. It is considered to be experimentally producible in high-energy 
heavy ion collisions at the Relativistic Heavy Ion Collider (RHIC) at the Brookhaven National 
Laboratory and the Large Hadron Collider (LHC) at the European Organization for Nuclear 
Research, providing great opportunities and challenges for developing realistic theories of the hot 
and dense QCD matter. One of the most distinctive characteristics of the QGP implied from the 
experimental data is that it behaves as a relativistic fluid with extremely small viscosity. This 
allows one to construct a hydrodynamic effective model of the high-energy nuclear collisions. 
Modern hydrodynamic studies incorporate viscosity and fluctuation for precision physics f2]. 

The net baryon number, so far, has been neglected in viscous hydrodynamic studies with an 
evolving flow. This could be due to the complexity of dissipative processes in causal relativistic 
formalisms and the boost invariant ansatz that fails to take account of the net baryon number 
conserved at forward rapidity. It should be noted that finite-density effects are generally more 
prominent near the beam rapidity since the net baryon number originates from the remnant of 
the shattered nuclei. Introducing the net baryon to the model would enable one to understand 
the matter-to-antimatter ratios of the particle spectra as well as the transport properties of 
finite-density media. The analyses of the finite-density matter via dissipative hydrodynamics 
would also be of great help in exploring the QCD phase diagram for the search of the critical 
point, given that first principle approaches have difficulties for calculating systems with a finite 
baryon chemical potential. As the beam energy scans are being performed and planned at RHIC, 
GSI, NICA and J-PARC, it is important to develop a finite-density model to analyze the data 
and to examine the validity of the model for those experiments. 



The net baryon distribution is quantified with baryon stopping defined by the mean rapidity 
loss {5y) = Vp — (y) [31I5- Here yp is the rapidity of the incoming projectile and 



is the average rapidity of the net baryon density. When {6y) = 0, the collision is completely 
transparent and no hot medium would be produced. On the other hand, {6y) = yp suggests that 
the nuclei comes to a complete halt at the collision, losing all the kinetic energy for the medium. 
Thus this quantify is related to the magnitude of the energy loss of the colliding nuclei that 
is available for the QGP production. Here I would like to focus on dissipative hydrodynamic 
effects on the net baryon distribution and the mean rapidity loss [5]. 

2. Dissipative hydrodynamic model 

The energy momentum tensor T^'^ and the net baryon current are the conserved quantities 
in the system of interest. Hydrodynamics introduces macroscopic flow which is defined here 
as the local energy flux. Then the tensor decomposition in terms of the flow yields 

T^^ = (eo + Po + n)n^n'^ - (Pq + ^)9^'' + = nBou^" + (2) 

where eo, Pq, nso, H, tt'^'^ and are the energy density, the hydrostatic pressure, the net 
baryon density, the bulk pressure, the shear stress tensor and the baryon dissipation current, 
respectively. The system is described by energy-momentum conservation d^T^^'^ = 0, baryon 
number conservation d^N^ = and the law of increasing entropy d^s^ > 0. By introducing 
second order corrections to the entropy current in terms of the dissipative quantities, the 
semi-positive definite condition yields causal hydrodynamic equations of motion [6l [7] : 



+ XWV^D^ + x'vvV^'D^ + xWy^'^uu'' + xUV'VuU^ + xWV'^^n, 



+ X^-nHV^^ + xUnV^^ + X^.nHI^n^' + X^nV^n, (4) 

+ x.nHV^'^^z-) + xl.^^'D^ + xL^'-'D^ + xl^^^''^ + xt^'^^ pW^^ 

+ xlvV^y^^Y' + ^-V'^^'^V-)! + xlvV^^Du''^ + xiv^^^y^. (5) 

Here D = u^d^ and = dfj, — Uf^D. The angular bracket on the indices denotes traceless 
symmetrization. T and fj,B are the temperature and the chemical potential, respectively. 
Ky and rj are the bulk viscosity, the baryon charge conductivity and the shear viscosity, which 
characterize scalar, vector and tensor dissipative processes in a medium at the linear order. (^uSe 



and Cn<5ns are the bulk-energy and the bulk-baryon cross viscosities and nyw is the baryon-heat 
cross conductivity, which should satisfy the Onsager reciprocal relations. It is note- worthy that 
they can be negative, rn, ry and are the relaxation times, which are required for causality 
and stability in relativistic systems, x are second order transport coefficients. 

One has to introduce the equation of state, the transport coefficients and the initial conditions 
to construct a model for the high-energy heavy ion collisions. The equation of state is obtained 
via Taylor expansion with the results from latest lattice QCD calculations IHUH]- The transport 
coefficients are obtained from the AdS/CFT approaches [lU^ 111), the non-equilibrium statistical 
operator method [12] and phenomenological analyses based on the dimensionality, the matter- 
antimatter symmetry and the semi-positive definiteness of the transport coefficient matrices [5j. 
It should be noted that the coefficients are introduced to investigate qualitative responses of 
the system. The initial conditions at tq = 1 fm are based on the color glass theory; the energy 
density is from the Nara Monte-Carlo version [13^ |T3] of Kharzeev-Levin-Nardi model \15\ [16] 
and the net baryon density from the distribution of the valence quarks \17\ I18j . The list of 
individual model input can be found in Tabled! where Xb quadratic baryon fluctuation, 

s is the entropy density and = [(9eo/5Po)s/nBo]^''^ sound velocity: 



Table 1. Input to dissipative hydrodynamic model 



Equation of state 


Description 


EoS at finite fis 


^o(eo, UBo) = Po(eo, 0) + {tilT^ /2)xf {eo, 0) 


Transport coefficients 


Description 


Shear viscosity 
Effective bulk viscosity* 
Baryon charge conductivity 
Baryon-heat conductivity 


rj — s/Aii 

Ceff = (5/2)[(l/3)-c2]r; 

i^v = cy(9/is/9ri_B)y"^/27r 

nvw = cvvF[»^Bo/(eo + fb)](577rKy)i/2 


Initial conditions 


Description 


Energy density 
Net baryon density 


eo(TO,??s) = {l/Se,rca)dET/Todr]s 

nBo{TQ,ris) = {I /S^^ca)dNg_g/ Todies 



*The contributions of the scalar cross coefficients are merged into the effective bulk viscosity. 

In the present analyses, I discuss dissipative hydrodynamic expansion in the longitudinal 
direction by integrating out the transverse geometry for the numerical analyses. It is note- 
worthy that the dependence of the net baryon yields on transverse geometry is shown to be very 
small experimentally in the analyses of centrality dependence |19j . 

3. Results 

Fig. [1] shows the net baryon rapidity distributions for the Au-Au collisions at RHIC and the Pb- 
Pb collisions at LHC for the most central 0-5 % events. The parameters for the conductivities 
are cy = 1 and cyw = for the moment. The contribution of hadron resonance up to 2 GeV 
is considered at the freeze-out temperature Tf = 0.16 GeV. The experimental data points are 
scaled results from the corresponding net proton distribution by the BRAHMS experiments 
[20]. One can see that the net baryon is carried to the forward rapidity region in hydrodynamic 
evolution. This is mainly due to the pressure gradient, to which the net baryon density is 
coupled via the flow. The effect is slightly moderate for the shear and bulk viscous case because 
the longitudinal pressure is effectively reduced by the viscosities. Baryon dissipation further 




Figure 1. The net baryon distributions of the initial condition and ideal, viscous and dissipative 
hydrodynamic results at RHIC (left) and at LHC (right) are shown as solid, dotted, dash-dotted and 
dashed curves, respectively. The experimental data points are from the BRAHMS Collaboration [20] . 



steepens the distribution as it is induced by the chemical gradient into the mid-rapidity region. 
For the current parameter settings, the hydrodynamic results for RHIC roughly agrees with the 
experimental data. The off-equilibrium corrections are more visible at RHIC than at LHC. 

The hydrodynamic evolution reduces the mean rapidity loss from ((5y)initiai = 2.67 to 
{Sy)hydm = 2.09,2.16 and 2.26 for the ideal, the viscous and the dissipative systems at RHIC. 
Here the dissipative results include effects of baryon diffusion in addition to those of shear and 
bulk viscosities. It shows that the transparency of the collision is sizably enhanced as a medium 
effect. This could explain the abrupt deviation of the RHIC data from the simple extrapolation of 
the AGS and SPS results (Fig. [2|) since the hot medium is created at the heavy ion collisions with 
sufficiently high energies. The kinetic energy losses after the hydrodynamic evolutions are also 
about 10-15% less than those of the initial conditions, suggesting that more energy is available 
for the QGP production at the time of collision than naively implied from the experimental 
data. At LHC, when (5y)initiai = 3.36, the rapidity loss is reduced to ((5y)hydro = 2.92 for the 
dissipative hydrodynamic case. It should be emphasized that the result is dependent on the 
choice of the initial conditions and the transport coefficients. 

Finally, the net baryon distributions in the presence of the baryon-heat cross-coupling term 
is shown in Fig. EJ The parameter for the cross coefficient is chosen as cwv = 5,0 and —5. The 
term induces Soret effect, a chemical diffusion process by the thermal gradient. The numerical 
results indicate that the effect is rather small. This is due to the fact that the cross conductivity 
has to stay small when the chemical potential is not large as Kvw{—fJ'B) = —i-^vwifJ-B) follows 
from the matter-antimatter symmetry condition V^^{—^b) = — ^'^(/^b) on Eq. 

4. Summary 

Non-equilibrium hydrodynamic analyses are extended for the quark-gluon systems with finite 
baryon density. The finite-density equation of state is constructed in the Taylor expansion 
method with state-of-art data of lattice QCD calculations and the initial conditions in the color 
glass-based pictures. The baryon stopping at RHIC and LHC show clear effects of hydrodynamic 
convection where both viscous and dissipative modifications are qualitatively visible at RHIC. 
The widening of the net baryon distribution suggests that the kinetic energy loss is effectively 
reduced during the hydrodynamic evolution. This implies that the energy available for the 
production of a hot medium is larger initially, and part of the energy is transferred back to the 




Figure 2. Mean rapidity losses at AGS [H], SPS 
[22] and RHIC [201 123] experiments and those in 
dissipative hydrodynamic calculations for RHIC 
and LHC. The dotted line is a linear extrapolation 
of the low energy results [50] . 



Figure 3. The net baryon distribution of the 
initial condition and dissipative hydrodynamic 
results at RHIC with the cross coefficient cvw = 
5, and —5 are shown in solid, dotted, dashed and 
dash-dotted curves, respectively. 



net baryon components at the later stage as the hydrodynamic flow accelerates them. Also the 
Soret effect would be small in the high-energy nuclear collisions. One would need to include 
transverse dynamics and realistic transport coefficients for more quantitative discussion. 
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